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Abstract. Let B H = {B H (t),t£ R+} be an (TV, d) -fractional Brownian sheet with index H = (Hi,.. .,H N ) £ (0, 1)^ 
defined by B H (t ) = (Bf (i) , . . . , Bf (t ) ) (t £ R+ ) , where 73f , . . . , Bf are independent copies of a real- valued fractional 
Brownian sheet Bq . We prove that if d < y\ H^ 1 , then the local times of B H are jointly continuous. This verifies 
a conjecture of Xiao and Zhang (Probab. Theory Related Fields 124 (2002)). 

We also establish sharp local and global Holder conditions for the local times of B H . These results are applied to 
study analytic and geometric properties of the sample paths of B H . 

Resume. Designons par B H — {B H (t),t £ R+ } le (TV, d)-drap Brownien fractionnaire de parametre H — (Hi, . . . , Hn) £ 
(0, 1)^ defini par B H (t) = (B"(t), B" (t)) (t G R+ ), oil Bf , . . . , Bf sont des copies independantes du drap Brow- 
nien fractionnaire a valeurs reelles Bq . Nous montrons que le temps local de B H est bicontinu lorsque d < X/^Li ■ 
Cela resout une conjecture de Xiao et Zhang (Probab. Theory Related Fields 124 (2002)). Nous obtenons aussi des 
resultats fins concernant la regularite Holderienne, locale et globale, du temps local. Ces resultats nous permettent 
d'etudier certaines proprietes analytiques et geometriques des trajectoires de B H . 
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1. Introduction 



For a given vector H = (Hi, . . . ,Hn) G (0, 1) N , a real-valued fractional Brownian sheet Bq = {Bq (t),t 6 



1+} with index H is a centered Gaussian random field with covariance function given by 



N 1 

E[BB(s)BB(t)] = J] M Hl + t] Hl fa t e \^), s,te: 



>N 



(1.1) 



It follows from (1.1) that Bq (i) = a.s. for every t 6 cM+ , where dM.^ denotes the boundary of ] 
We will make use of the following stochastic integral representation of Bq (cf. [2]): 



B (*) - K H 



tN N 

Y[g He (t e , Se )W(ds), 

■ 00 i=i 



(1.2) 



where W = {W(s),s G R w } is a standard real- valued Brownian sheet and where, for every 

g Bt (U, 8i ) = ({t e - s e ) + ) H '- 1/2 - ((- s ,) + ) ff - 1/2 . 
In the above, a+ = max{a, 0} for all a £ M and k h is the normalization constant given by 



1,- 



,N, 



.£=1 



ds. 



Note that if H^ = 1/2 for some £q, then we assume that g H (t( a ,se ) = i[o,t £ ] (s^ ); where l[o.t £ ] is the 
indicator of the interval [0,^ o ]. 

Let , ...,Bf bed independent copies of Bff . Then the Gaussian random field B H = {B H (t): t £ R^} 
with values in W l defined by 



B H (t) = (B?(t),...,B?(t)), Vie 



(1.3) 



is called an (TV, d) -fractional Brownian sheet with index H = (Hi, . . . , Hn). 

Note that if N = 1, then is a fractional Brownian motion in E d with Hurst index Hi € (0, 1); if N > 1 
and Hi = ■ ■ ■ = Hm = 1/2, then is the (iV, d)-Brownian sheet. However, when Hi, ... , Hn are not the 
same, _B ff is anisotropic and has the following operator-self-similarity (this can be verified easily using (1.1)): 
For any N x N diagonal matrix A = (ay) with an = a% > for all 1 < i < N and ay = if i ^ j , we have 



JV 



{B H (At),te^}=l'[[af i B H (t),t&: 



(1.4) 



where X = Y means that the two processes have the same finite dimensional distributions. These features 
of B H make it a possible model for bone structure [8] and aquifer structure in hydrology [4]. 

Many authors have studied various properties of fractional Brownian sheets. See, for example, [3, 11, 
21, 24, 29, 33] and the references therein for further information. This paper is concerned with regularity 
of the local times of an (N, d)-fractional Brownian sheet B H . After having proved that a necessary and 
sufficient condition for the existence of L 2 (P x Ad) local times of B H is d < TT' Xiao and Zhang [33] 

give a sufficient condition for the joint continuity of the local times. However, their sufficient condition is 
not sharp and they have conjectured that B H has jointly continuous local times whenever the condition 
d < Yle=i 777 * s sa ti sn °d. The main objective of this paper is to verify this conjecture; see Theorem 3.1. The 
new ingredients for proving this result is the property of sectorial local nondeterminism of Bq established 
in [29] (see Lemma 3.2) and a similar result for the fractional Liouville sheet proved in Section 2. The results 
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and techniques developed in this paper are applicable to more general anisotropic Gaussian random fields 
with the property of sectorial local nondeterminism; see [32] for further development. 

The rest of this paper is organized as follows. In Section 2, we prove some basic results on the fractional 
Liouville sheets that will be useful to our arguments. In Section 3, we prove that the sufficient condition for 
the existence of L 2 (P x Xd) local times of B H in [33] actually implies the joint continuity of the local times. 
This verifies their conjecture in Remark 4.11. Section 4 is on the local and uniform Holder conditions for the 
local times and their implications to sample path properties of B H . In particular, we derive some results on 
the Hausdorff measure of the level sets and on the Chung-type law of the iterated logarithm for the sample 
function B H (t). The latter improves Theorem 3 of [3]. 

We end the Introduction with some notation. Throughout this paper, the underlying parameter space is 
'R N or M.+ = [0, oo)^. A parameter t £ M. N is written as t = (ti, . . . , t^), or as (c), if t\ = ■ ■ ■ = tjy = c. For any 
s,t € l w such that Sj <tj (j = 1, . . . , N), we define the closed interval (or rectangle) [s,t] = Y\J =1 [sj ,tj]. We 
will let A denote the class of all closed intervals / C (0, oo) N . We always write A m for Lebesgue's measure on 
R m , and use (•, •) and | • | to denote the ordinary scalar product and the Euclidean norm in R m respectively, 
no matter the value of the integer m. 

An unspecified positive and finite constant will be denoted by c, which may not be the same in each 
occurrence. More specific constants in Section i are numbered as c i ls c 4 2 , 

2. Fractional Liouville sheet 

One of the main obstacles in studying the local times and other properties of fractional Brownian sheets is 
their complicated dependence structure. Unlike the Brownian sheet or fractional Brownian motion, fractional 
Brownian sheets have neither the property of independent increments nor the local nondeterminism. 

To be more specific, we recall that fractional Brownian motion Z a = {Z a (t),t £ 1™} (0 < a < 1) in K 
has the following property of strong local nondeterminism proved by Pitt [25]: For every interval I C M. N , 
there exist positive constants c 2 a and r such that for alH £ / and all < r < min{|£|, ro}, 

Var(Z Q (f) | Z a (s): s e I ,r < \s - t\ < r ) > c 21 r 2a . (2.1) 

This property has played important roles in studying the local times and many other properties of Z a ; see 
[31] and the references therein for more information. On the other hand, it is known that the Brownian 
sheet W = {W(t),t G R^} does not have the property of local nondeterminism. In order to see this, we 
consider the Brownian sheet with N = 2 and I = [0, l] 2 . For any constant e £ (0, 1), let T C / be an interval 
with side-length e. Let t denote the upper-right vertex of T and let s^s^s 3 be other vertices of T. For 
example, t = (1,1), s 1 = (1 - e, l),s 2 = (1,1 - e) and s 3 = (1 - e, 1 - e). Then \t - s j \ > e for j = 1,2,3. 
Considering the increment of W over the square T, we see that Va.r (W(t)\W(s 1 ), W(s 2 ), W(s 3 )) < e 2 . Hence 
the Brownian sheet W does not satisfy (2.1) (this also proves that the fractional Brownian sheet B$ is not 
locally nondctcrministic). This is the main reason why, in most literature, the methods for studying various 
properties of the Brownian sheet are different from those for fractional Brownian motion. The property of 
independent increments of W has been crucial in studying the local times and self-intersection local times 
of W; see [12, 27] and [22]. In solving an open problem in [22], Khoshnevisan and Xiao [19] showed that W 
satisfies a type of sectorial local nondeterminism and applied this property to study geometric properties of 
the Brownian sheet by using methods that are reminiscent to those for fractional Brownian motion; see [18] 
for further applications of the sectorial local nondeterminism. Recently, Wu and Xiao [29] have extended 
several results in [18, 19] to fractional Brownian sheets. 

In this paper we continue the above line of research and study the regularity of the local times of fractional 
Brownian sheets. To overcome the difficulty due to the lack of local nondeterminism of B H , we will not only 
make use of the sectorial local nondeterministic property of B H established in [29] (see Lemma 3.2), but 
also the analogous properties of the so-called fractional Liouville sheet. 
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Given any vector a = (ai, . . . , a/v) £ (0, oo) , the centered Gaussian random field Xq = {Xq (t),t £ M^} 
denned by 

N 

*o(*)=/ l\(t e - S t) a *-V 2 w(ds), ten%, (2.2) 

J[o,t] e=1 

is called a fractional Liouville sheet with parameter a. It is easy to see that, when ai, . . . , ajv are not the 
same, Xq — {Xq (i): t £ R^} is an anisotropic Gaussian field which has the same operator self-similarity as 
in (1.4). 

For the purpose of this paper, we will only be interested in the case a = H £ (0, 1)^. It follows from (1.2) 
that for every t £ R+ , 

N 

B Q H {t)=n- H 1 X H {t) + K H 1 Jg (t £ ,s e )W(ds), (2.3) 

J(.-oo,t]\[0,t] t _. 



and the two processes on the right-hand side of (2.3) are independent. We will show that in studying the 
regularity properties of the local times of B H , the fractional Liouville sheet X^ plays a crucial role and the 
second process in (2.3) can be neglected. More precisely, we will make use of the following property: For all 
integers n > 2, t 1 , . . . , t n € R^ and «i,...,u„£R, we have 

Var (j2 %W)J > «a Var ^ UjX^)j . (2.4) 

Here and in the sequel, Var(^) denotes the variance of the random variable £. 

Next we use an argument in [3] to provide a useful decomposition for Xq (t). Let e > be fixed. For every 
t € [e,oo) , we decompose the rectangle [0,t] into the following disjoint union of sub-rectangles: 

N 

[0,t] = [0,e} N u\jR e (t)UA(e,t), (2.5) 
t=\ 

where Ri(t) := Rg(e,t) ={r£ [0,^^: < r; < e if i ^ £, £ < ri < t^} and A(e,t) can be written as a union 
of 2^ — N — 1 sub-rectangles of [0,t]. Denote the integrand in (2.2) by g(t,r). It follows from (2.5) that for 
every t £ [e, oo) N , 

X H (t)= f g(t,r)W(dr) 

J[0,e] N 

+ Y] I g(t,r)W(dr)+ [ g(t,r)W(dr) 

e=1 JR e (t) JA(e,t) 
N 

:=X(e 7 t) + J2Ye(t) + Z(e,t). (2.6) 
l=i 

Since {X(e,t),t£ [£,00)*}, {Y e {t),t£ [£,00)^} (1 < £ < N) and {Z{e,t),t£ [£,00)^} are defined by the 
stochastic integrals w.r.t. W over disjoint sets, they are independent Gaussian random fields. 

The following lemma shows that every process Yi(t) has the property of strong local nondeterminism 
along the £th direction. It will be essential to our proofs. 

Lemma 2.1. Let i £ {1, 2, . . . , N} and let I = [a, b] £ A be a fixed interval. For any integer n > 2, t , . . . , t n £ 
[a, b] such that 

t)<ti <•••<*?, 
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the following inequality for the conditional variance holds: 

Vfa{Yt(t n )\Y t (lt): 1 < j < n - 1) > c 2 2 \t n t - t n ~ x \ 2H \ (2.7) 
where c 2 2 > is a constant depending on e, I and H only. 

Proof. Working in the Hilbert space setting, the conditional variance in (2.7) is the square of the L 2 (P)- 
distance of Yi{t n ) from the subspace generated by Yg(P) (l<j<n — 1). Hence it is sufficient to show that 
there exists a constant c 2 2 such that for all aj G R (j = 1, . . . , n — 1), 

E^(n-E a ^( <3 )) -tr x \ 2He - (2-8) 

However, by splitting Ri(t n ) into two disjoint parts and using the independence, we derive that 

2 \ 2 
g(t n ,r)W(dr) 



Jo Jt" 1 Jo 



k=l 

>c,,j*?-*rT* £ - (2-9) 

This proves (2.8) and hence Lemma 2.1. □ 

The following lemma relates the fractional Brownian sheet Bq to the independent Gaussian random 
fields Yj> (£ = 1,...,N). 

Lemma 2.2. Let I = [a, b] G A be a fixed interval. For every integer n > 2, t , . . . , t n G [a, b] and Ui, . . . , u n € 

R, we have 

(n \ N / n \ 

J2 UjBgtf) > Kff 2 Var ^ % W) . (2.10) 
i=i / <=i \j=i / 

Moreover, for every k G {1, . . . , N} and positive numbers pi, . . . ,pk > 1 satisfying "Y^l^iP^ 1 =1, we have 

1 k c n 

[detCov(B? (t 1 ), . . . ,£<f (^))]V2 " 2 [dctCov(r,(ti), . . ., (*"))] 1/(2pf) ' (2 ' U) 

where detCov(Zi, . . . , Z n ) denotes the determinant of the covariance matrix of the Gaussian random vector 
(Zi, . . . , Z n ). 

Proof. The inequality (2.10) follows directly from (2.4), (2.6) and the independence of Ye (£= 1,...,N). 
To prove (2.11), we note that for any positive definite n x n matrix r, 

[det(r)]V» / 1 



/• [det(i 



exp --iTi <k = l. (2.12) 



t)"/ 2 

It follows from (2.10), (2.12) and the generalized Holder inequality (see, e.g., [15], p. 140) that 
1 



[detCov(B?(ti),.„,B?(t»))]Va 
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(2 



< 



dp72 X« exp ~5 Var (£ u ^^J dM i'" du « 

N / n \ 

c E Var E^W) 



(27t)™/ 2 



< 



< 



(271) 



n 



cxp 



\j=i > 



dux ■ • • du n 



dui ■ ■ ■ du r 



^ [detCov(r £ (ii), . . . ,y*(* n ))] 1/(2w) 

This yields (2.11) and the lemma is proved. 



(2.13) 
□ 



3. Joint continuity of the local times 

We start by briefly recalling some aspects of the theory of local times. For excellent surveys on local times 
of random and deterministic vector fields, we refer to [13] and [10]. 

Let X(i) be a Borel vector field on R w with values in R d . For any Borel set T C R , the occupation 
measure of X on T is defined as the following measure on R d : 

M») = Ajv{t S T: *(<)€•}. 

If /i T is absolutely continuous with respect to A^, we say that X(t) has local times on T, and define its 
local times, L(*,T), as the Radon-Nikodym derivative of n T with respect to A<j, i.e., 



L(x,T) 



dA d 



(x), VxeR d . 



In the above, x is the so-called space variable, and T is the time variable of the local times. Sometimes, 
we write L(x,t) in place of L(x, [0,t]). Note that if X has local times on T then for every Borel set SCT, 
I/(x, «S) also exists. 

By standard martingale and monotone class arguments, one can deduce that the local times have a 
version, still denoted by L(x,T), such that it is a kernel in the following sense: 

(i) For each fixed S € 15(7), where 25 (T) is the family of Borel subsets of T, the function x i— > i(x, S 1 ) is 
Borel measurable in x S R d . 

(ii) For every x £ R d , L(x, •) is Borel measure on 25 (T). 

Moreover, L(x,T) satisfies the following occupation density formula: For every Borel set T CR N , and for 
every measurable function / : R d — > R + , 



/(*(*)) dt= / /(x)L 



(x, T) dx. 



(3.1) 



See Theorems 6.3 and 6.4 in [13]. 

Suppose we fix a rectangle I = Il^Li[ a i> a i + hi] m Then, whenever we can choose a version of the 
local time, still denoted by L(x,Y\f =1 [ai : ai +U]), such that it is a continuous function of (x, ti, . . . , tjv) 
€ R d x n i=1 [0, /ij], X is said to have a jointly continuous local time on /. When a local time is jointly 
continuous, L(x, •) can be extended to be a finite Borel measure supported on the level set 



X- 1 (x)n/ = {te/: X(t)=x}; 



(3.2) 
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see [1] for details. In other words, local times often act as a natural measure on the level sets of X. As such, 
they arc useful in studying the various fractal properties of level sets and inverse images of the vector field 
X. In this regard, we refer to [6, 12, 27] and [30]. 

Berman [5, 6, 7] developed Fourier analytic methods for studying the existence and regularity of the 
local times of Gaussian processes. His methods were extended by Pitt [25] and Geman and Horowitz [13] to 
Gaussian random fields. Let X = {X(t),t G 1H N } be a Gaussian random field with values in M. d . It follows 
from (25.5) and (25.7) in [13] (sec also [14, 25]) that for all x, y G R d , T G A and all integers n > 1, 



E[L(x,Ty i ] = (2n)- nd [ [ expf-iV^,: 

J T nJ R nd V J 



x Ecxp^i^(w J ,A(i J ))J dudt (3.3) 



and for all even integers n > 2 , 

E[(L(x,T)-L(y,T)) n ] = {2n)- nd J J Jjje" 1 ^ - e" 1 ^] 



Eexp^i^^'.X^J dudt, (3.4) 



where u = (u 1 , . . . ,u n ),t = (t 1 , . . . ,t n ), and each £M. d ,P eTc (0,00)^. In the coordinate notation we 
then write u J — (u\, . . . These identities are also very useful for studying the local times of infinitely 
divisible random fields as well; see [10, 12] and [20]. 

Xiao and Zhang [33] have proved that if d < J2i=i then for all intervals I G A, B H has local times 
{L{x,I),x G K d } on / and L(-,I) G i 2 (P x Ad). In the following, we prove that under the same condition, 
the local time has a version that is jointly continuous in both space and time variables. 

Theorem 3.1. Let B H = {B H (t),t G R^} be a fractional Brownian sheet in M. d with index H = 
{H\, . . . , Hn) G (0, 1) N . If d < X^=i Jf~i then for all intervals I & A, B H has a jointly continuous local 
time on I almost surely. 



To prove Theorem 3.1 we will, similar to [12, 30, 33], first use the Fourier analytic arguments to derive 
estimates on the moments of the local times (see Lemmas 3.7 and 3.10) and then apply a multiparameter 
version of Kolmogorov continuity theorem (cf. [17]). The new ingredients in this paper are the "sectorial 
local nondeterministic" properties of fractional Brownian sheets proved in [29] and the results on fractional 
Liouville sheets proved in Section 2. 

We will also make use of the following lemmas. Among them, Lemma 3.2 is proved in [29] and Lemma 
3.3 is essentially due to Cuzick and DuPrcez [9] (see also [19]). 

Lemma 3.2. Let Bq = {Bq (t), t G M^} be a fractional Brownian sheet in K with index H = {Hi, . . . , Hm) G 
(0, 1) N . Then for every e > 0, there is a constant c 3 1 > such that for all integers n > 2, t 1 , . . . , t n G [e, 00) , 

N 

Var«(r) I B^ (a*), j + n)> c 3>1 X>in{|t? " < j<n- 1}, (3.5) 

1=1 



where £? = for I = 1, . . . , N. 
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Lemma 3.3. Let Z\, . . ., Z n be mean zero Gaussian variables which are linearly independent, then for any 
nonnegative Borel function g : K — > K. + , 



dvi ■ ■ ■ dv„ 



(detCov(Z 1 ,...,Z,M l l* i_„ 9 V»i 
where o\ = Vai(Zi\Z2, ■ ■ ■ , Z n ) is the conditional variance of Z\ given Z2, . . . ,Z n . 

The following technical lemma is essential in establishing the moment estimates for the local times L(x,T). 
Since it may be of independent interest, we state it in a more general form than is needed in this paper. 

Lemma 3.4. For any q £ [0, J^&Li ^i" 1 )* let t G {1, . . . , N} be the integer such that 

with the convention that Y2i=i jij : ~ 0- Then there exists a positive constant S T < 1 depending on 
(Hi, . . . , Hn) only such that for every S € (0, 5 T ), we can find t real numbers pg>l (1 < I < t) satisfy- 
ing the following properties: 



E- = l, — <L V/=l,...,r, (3.7) 
j^Pt Pi 

and 

{1 . S) ±^l £HTq + T .±f. ,3. 8) 

Furthermore, if we denote a T := 2l=i jf ~ 9 > 0> i/ien /or anj/ positive number p £ (0, ^jr)> i/iere exists an 
index £q £ {1, . . . ,r} smc/i £/iai 

^ + 2H loP <l. (3.9) 

Wo 

Remark 3.5. It is important to note that the choice of the numbers pi>l (1 <£<t) depends on 5. 
Moreover, it follows from the proof below that, except for the case of t = 2, we can always take 5 T = 1. 

Proof of Lemma 3.4. First we prove (3.7) and (3.8). If (3.6) holds for r = 1, then for all < S < Si := 1, 
we can take Pi = l and both (3.7) and (3.8) hold automatically. 

We now prove the cases of r > 2 by induction. Our proof provides a general procedure for constructing a 
sequence {pe, 1 < £ < t} of real numbers pi>\ satisfying (3.7) and (3.8) (there are many possible choices). 

Assume that (3.6) holds for r = 2. We distinguish two cases: (i) Hi = H2 and (ii) Hi ^ H^- In the first 
case, we have Hi^ 1 < q < 2Hi~ 1 . We choose 77 > such that 

(2-H iq )H iq 

0<V< -r ; 

Hiq - 1 

(if Hiq = 1, then n > can be arbitrarily chosen) and define 

— = and — = 1 . 

pi Hiq + r) p 2 Pi 
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Then a few lines of calculation verify that p\ and P2 satisfy (3.7) and (3.8) for all 8 £ (0, 1). 

To consider the case (ii) we may and will assume, without loss of generality, that Hi < Hi- Since q < 
H^ + H2 , there exists 82 > such that for all S £ (0, 62), 

H 1 H 2 q(H 2 -H 1 +SH 1 )< (fla-lfi)(ifa+-ffi - SH ± ). (3.10) 
For each fixed 8 £ (0, 62), we define 

1 1 1 8 H 2 ad 1 -! 1 

pi 1 — 8 H\q 1—8 H2 — Hi P2 Pi 

Then (3.7) follows from (3.6) and (3.10), and the equality sign in (3.8) holds. 

Now we assume that the properties (3.7) and (3.8) hold for t = n £ {2, . . . , TV — 1} and consider the case 
of t = n + 1 . Then we have 

n . . n+1 1 

1=2 1=2 

Let 8 £ (0,1) be fixed and we choose 8' £ (0,8 A 8 n ). Then it follows from (3.11) and the induction 
hypothesis that there exist n constants p\>l (£ = 2, . . . ,n + 1) such that 

£1 = 1, (g ~ 1/gl) ^ <l, W = 2,...,n + 1, (3.12) 

and 

^H*(q-l/HJ ( 1\ n ^H n+1 

1=2 ^ l v L/ 1=2 1 

To define the constants pi, . . . ,p n +\ with the desired properties, we choose a constant rj > small so that 



Heq( l 1 



o'e \ Hiq 
and 

(l-a)(l + (ffig77)/(ffig-l)) 

i-<y 



t? <1, W = 2,...,n + 1, (3.14) 



<1. (3.15) 



This is possible because of (3.12). 
Now wc define pi (1 < I < n + 1) by 

i = l(l_J-+J, W = 2,...,n+1, (3.16) 

and 

- = -=--»,. (3.17) 

Pi Hiq 

It follows from this definition and (3.14) that 

n+l 



n+1 „ 

£ — = 1 and — <1, W=l,2,...,n+1. (3.18) 



e=i 
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That is, (3.7) holds for r — n + 1. On the other hand, by some elementary calculation and (3.15) we can 
verify that 

n+1 „ n+1 „ 

(l-^7< H n+iq +(n + l)-J2^. (3.19) 

That is, (3.8) also holds for r = n + 1. Hence the proof of (3.7) and (3.8) is completed. 

Finally we prove (3.9). By (3.7), for every £ £ {1, r}, 3ee € (0, 1) such that ^ = l- £ ^. Hence, 

^ He ^ He ^ p? ^ Hp 

i=i t=i i=i ^ t=i 



Hence there exists ^6 {1, • ■ • , t} such that 6( > g 'f T ■ Note that for every positive number p € (0, y l ), we 
have 2i^ oj o < < £4. Therefore 

^+2i^ p = l-^ + 2tf^<l, (3.21) 

Pio 

which completes the proof of (3.9). □ 

The following inequalities (3.22) and (3.23) with a = arc well known; see, e.g., [12]. The case a > 
makes it possible for us to apply Lemma 3.4 for proving Lemmas 3.7 and 3.10. 

Lemma 3.6. For all integers n>l, positive numbers a, r, < bj < 1 and an arbitrary so G [0,a/2], 

/ ~[[( Sj - Sj .i)-^d Sl ... ds n <cl a (r*f /n) ZU b i-\"-i:U b i, (3.22) 

J a<si < ■ ■ -<s n <a+r j — \ 

where c 3 2 > is a constant depending on a and bj J s only. In particular, if bj = a for all j = 1, . . . , n, then 
/ n(s J -s J - 1 )- a dsi---ds n <cl 2 (n\) a ~ 1 r n ( 1 -( 1 - 1 / n ^. (3.23) 

J a<si<---<s n <a-\-r j = \ 

Proof. For simplicity, we only give the proof of (3.23) here. The proof of (3.22) is almost identical, and thus 
omitted. By integrating the integral in (3.23) in the order of ds„,ds„-i, . . . ,dsi, and by using a change of 
variable in each step to construct Beta functions, we derive 

« n 

/ ~[( s 3 ' -Sj-i)~"dsi--- ds n 

In— 2 ra+r 



i r(2-a)[r(i-a)] 1 



(a + r- ai )(n-i)(i-a) ( Sl _ 8Q ya dsi (3 24) 



1-a r(l + (n-l)(l-a)) J a 
The inequality (3.23) follows from (3.24) and the Stirling's formula. □ 

In the rest of this section, we assume that d < Yl(.=i TT e aim I G A is a fixed interval. For convenience, we 
further assume in the rest of this paper that 

0<H 1 <---<H N <1. (3.25) 

We proceed to establish the moment estimates for the local times L(x,T) which will be useful for proving 
the joint continuity of local times. 
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Lemma 3.7. Let B H = {B H (t),t G R^} be a fractional Brownian sheet in M. d with index H = (Hi, . . . , Hm) 
If for some integer r € {1, . . . , N} we have 



r-l t 1 

f-*> Hi~ ^ Ho 



(3.26) 



f ften 4/iere exists a positive and finite constant c 3 3 , depending on N, d, H and I only, such that for all 
intervals T = [a, a + (r)] C 7 wit/i edge-length r E (0, 1) , aZZ x S R d and a?Z integers n > 1, 



E[L(x,T) n ] < c n i3 (n\) N -^r n ^ , 
where /3 t =N-t- H r d + ££ =1 ff T /F<> 



(3.27) 



Remark 3.8. As we mentioned earlier, the local time L(x,») may be extended as a random Borel measure 
supported on the level set F x = {t 6 (0, oo) : B (t) — x} . Hence the moment estimate (3.27) contains a lot 
of information about the fractal properties of r x . By Theorem 5 of [3], the Hausdorff dimension of the level 
set is given by 



dim H r x = min 



H, 



N-k-H k d + J2jfA<k<N\, 



e=i 



(3.28) 



and the minimum is achieved by j3 T = N — r — H T d+ Yle=i Ht/Hi, where t satisfies (3.26). It is important 
to note that (3.27) is sharp and can be applied to strengthen the Hausdorff dimension result (3.28). We 
believe that the function <fii(r) = r 3x (loglog l/r) N ~@ r is an exact Hausdorff measure function for r x , and 
we will give a proof for the lower bound of the ipi- Hausdorff measure of T x in Section 4- However, since the 
upper bound part relies on different methods, we will have to deal with it elsewhere. 

Proof of Lemma 3.7. For later use, we will start with an arbitrary closed interval T = YieLi i a e, a-e + r?} C I. 
It follows from (3.3) and the fact that B^ , . . . , B^ are independent copies of Bq that for all integers n > 1, 

E[L(x,T)"]<(27t)-" d ^n| j („ ex P -5 Var ^X(*'')j dlljjdt, (3.29) 

where Uk = (ui, ■ ■ ■ , u%) € K™. Fix k = l,...,d and denote the inner integral in (3.29) by 3k- Then by Lemma 
2.2, we have 



2k < 



cxp 



< 



i«H 3 EVar(£y fc y^) 
e=i \j=i j 

f exp -^ H 2 E Var (E^^') 



dU k 



dU k . 



(3.30) 



Since (3.26) holds, we apply Lemma 3.4 with S = n^ 1 and q = d to obtain t positive numbers pi, . . . ,p T > 1 
satisfying (3.7) and (3.8). 

Applying the generalized Holder inequality ([15], p. 140) to the last integral in (3.30), wc derive that 



a ' £ S{/„ 



-y K H 2var ( E u fe y ^* 
JidetCov^t 1 ),...^))} 



3=1 / 

n\\-i-l/(2pe) 



I/Pt 



dU k 



(3.31) 



738 



A. Ayache, D. Wu and Y. Xiao 



where the last equality follows from (2.12). Hence it follows from (3.29) and (3.31) that 

E[L(x,T) n ]<cl 5 f f[[detCov(Y,(t 1 ),...,Y i (t n W d/i2Pi) ^. (3.32) 
Jt™ i=1 

To evaluate the integral in (3.32). we will first integrate [dtj ■ ■ ■ dt™] for I = 1, . . . , r. To this end, we will 

make use of the following fact about multivariate normal distributions: For any Gaussian random vector 
(Zi, . . . , Z n ), 

n 

detCov(Zi, . . . , Z n ) = Var(Zi) JJ Var(^ |Zi, . . . , Z^). (3.33) 

3=2 

By the above fact and Lemma 2.1, we can derive that for every t £ {1, . . . , r} and for all t 1 ,. . . , t n £ T = 
Y\f =1 [a e ,a e + r e ] satisfying 

ai < t; e{1) < tj l{2) <■< t; ein) <a e + r £ (3.34) 
for some permutation tt£ of {1, . . . , AT}, we have 

n 

detCov^i 1 ), . . .MH) > < 6 U(t? (j) ~ tf^f 11 ', (3.35) 

i=i 

where tj'^ :— e (recall the decomposition (2.6)). We have chosen e < \ min{a^ , 1 < £ < N} so that Lemma 
3.6 is applicable. 

It follows from (3.34) and (3.35) that 



f [detCovO^ 1 ),...,!^"))] 



n ^- d/i2Pe) dt) ■ ■ ■ dtl 



< y c " / t — n — p— dt\ ■ ■ ■ dtf 

< C™ ( n \} H e d /Pe r "( 1 -( 1 - 1 / n ) H e d /P<!) _ (3.36) 

In the above, the last inequality follows from (3.23). 

Combining (3.32), (3.36) and continuing to integrate \dt\ ■ ■ ■ dt™) for i = r + 1, . . . ,N, we obtain 

T N 

E[L(x,r) n ]<c" 8 (n!)^=i^ ti/Pf TJ^ (1 ~ (1_1/n)fffd/w) - If r "' ( 3 - 37 ) 

1 = 1 l=T + l 

Now we consider the special case when T = [a, a + (r)], i.e. r-y = ■ ■■ = rjv = r. Equations (3.37) and (3.8) 
with 6 = n^ 1 and q — d together yield 

E[L{x,T) n ] < c^ 9 (nl)^=i Hid/pi r n{N - {1 - n ' 1) ^=i Hed/pe) 

<c n 3iQ {n\) N -^r n ^. (3.38) 
This proves (3.27). □ 

Remark 3. 9. In the proof of Lemma 3. 7, if we apply the generalized Holder inequality to the first integral 
in (3.30) with N positive numbers pi, ■ ■■ ,Pn defined by 

N H 

W = Ep 1=1,.. .,N, 
i=i 1 
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then the above proof leading to (3.37) shows that the following inequality 
E[L(x, T) n ] < <£ u {n\) Nv X N {T) n ^-^ 



(3.39) 



holds for every interval T C i\ where v = d/ 1 ) G (0, 1). We will apply this inequality in the proof 
of Theorem 3.1. 

Lemma 3.7 implies that for all n > 1, L(x,T) G L n (W d ) a.s. (see [13], p. 42). Our next lemma estimates 
the moments of the increments of L(x,T) in x. 

Lemma 3.10. Assume (3.26) holds for some t G {1, N} . Then there exists a constant c 3 12 , depending 
on N, d, H and I only, such that for all hypercubes T = [a, a + (r)} <Z I , x, y G M. d with \x — y\ < 1, all even 
integers n > 1 and 7 G (0, 1 A §^), 



E[(L(x,T) - L(y,T))"] < c n 3 12 {n\) N -^ +{1+H ^\x - y \^ r n ^- H ^\ 
Proof. Let 7 G (0, 1 A be a constant. Note that by the elementary inequalities 

| e i«_l| <2 1 -T| U |T foraUuGM 
and |u + 1>| 7 < |u| 7 + |t>| 7 , we see that for all u 1 , . . . , u n , x, y G M. d , 



■it, it 



(3.40) 



(3.41) 



(3.42) 



3=1 



where the summation is taken over all the sequences (ki, . . . , k n ) G {1, . . . , d} r 
It follows from (3.4) and (3.42) that for every even integer n > 2, 

E[(L(x,T)-L(y,T)) n ] 
< (2n)- nd 2 {1 -"< )n \x ~ y\ n ~< 

C C n / n \ 

x E / / I K m J 7E H -iE^' sff ^)) d « d * 
<0-yrE'/ dZ 

1 / - 

--Var E^'.-BV)) 



n 



du 



(3.43) 



where the last inequality follows from the generalized Holder inequality. 

Now we fix a vector k = (ki, fe, • • . , k n ) G {1, . . . , d} n and n points t ,. . . ,t n G T such that t\, . . . , i" are 
all distinct for every 1 < I < N (the set of such points has full (nA^)-dimensional Lebesgue measure). Let 
M = M(k,t,j) be defined by 



M =ML 



,m i™7 



exp 



1/n 



du 



(3.44) 



Note that Bf (1 < i < N) are independent copies of Bq . By Lemma 3.2, the random variables Bf{V) 
(1 < I < N, 1 < j < n) are linearly independent. Hence Lemma 3.3 gives 



1 m 1 Tvy 

l«L ex P 



di/ 
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(27t )(nd-l)/2 f f v \ nl 



" [detCovOB^ 1 ),..., 5<f (t"))] d / 2 J R \a, 
c" i4 (n!)T i 

" [detCov(£*(t"), . . . , B* (i«))] d/2 (3 ' 45) 

where o - ^ is the conditional variance of (t m ) given B^ (tP) (i^k m or i — k m but j ^ m), and the last 
inequality follows from Stirling's formula. 
Combining (3.44) and (3.45) we obtain 



M< 



c n (n^V n 1 

n ±- <3«) 



[detCov(S "(<l) BS(P))\" n ii, «l 



For 6 = l/n and q = d, let (£ = 1,...,t) be the constants as in Lemma 3.4. Observe that, since 
7 £ (0, § 3L ), there exists an £ G {!,-.., such that 



H £a d 



2H iol < 1. (3.47) 

Wo 

It follows from (3.46) and Lemma 2.2 that 

*<<^ w-n [de.covm*-), 1 .,^))]^-.' fig- < 3 - 48 ' 

The second product in (3.48) will be treated as a "perturbation" factor and will be shown to be small 
when integrated. For this purpose, we use again the independence of the coordinate processes of B H and 
Lemma 3.2 to derive 

a 2 m =V & r(B?Jt m )\B£jti),j^m) 

N 

> c 2 17 £min{|t? - 4\ 2H ': j ± m}. (3.49) 



For any n points t 1 , . . . , t n 6 T, let tt\ , . . . , ttn be N permutations of {1, 2, . . . , n} such that for every 1 < I < N, 
<e/ 2) <---<e/ (n) . (3.50) 
Then, by (3.49) and (3.50) we have 



n 1 n 



11 — ; 

m— 1 



1 



= l C 3,isH t £ c £ i A 1% % 7J 

- 1 

- C 3™8 II , W<0 (m) ,^ (m-l).,".ff, o7 ' ( 3 - 51 ) 

for some (g] o , . . . , g™ ) € {0, 1, 2} n satisfying J2m=i 1e^ = n anc ^ 1i = ^- That is, we will only need to consider 
the contribution of a m in the ^oth direction. 
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So far we have obtained all the ingredients for bounding the integral in (3.43) and the rest of the proof 
is quite similar to the proof of Lemma 3.7. It follows from (3.48) and (3.51) that 



1 



L JL [detCov(y £ (t 1 ), . . .,Y e (t n ))] d /( 2 P*) 

n 

m=l \i-e % ) 

To evaluate the above integral, we will first integrate [dt], ■ ■ ■ dt™] for every I = 1, . . . , r. Let us first consider 
I = £q. By using Lemma 2.1, (3.33), (3.22) and, thanks to (3.47) and the nature of q^, we see that 

Jla to ,a eo+ r to] n [detCov(y, (tl),...,y, (t™))] d /(^o) X ^ ( ^ (m)_^ (m-l) ) , roH<o7 d 4 ' ' ' dt "o ( 3 ' 53 ) 

<y c » [ 

— L^l 3,20 / « 4 1 

n 

x n (C (m) - ^ (m - i) )- ( ^° dMo+9r °^ o7) < • ■ • d^ o 

m— 1 

<c™ ( n !) ff fo d /Pf + ff fo'r r ™ [1 ~ (1_1/,l)fffod/p,!o " ff,!o71 . (3.54) 

In the above, tj e °^ = e as in the proof of Lemma 3.7 and the last inequality follows from (3.22). 
Meanwhile, recall that, for every £ ^ £q (£ S {1, . . . , t}), we have shown in (3.36) that 



/ [detCov^^ 1 ),...,^^ 1 ))]- 

J \a/,a/ +r»l" 



< c ™ 7 („!)^ d /PV™ (1 " (1 " 1/n) ^ d/pf) . (3.55) 



Finally, we proceed to integrate [dt] • • • dt"] for I = r + 1, . . . , N. It follows from the above that 
/ M(fc,t, 7 )dt < c " 22 (n!)£I=i fffd/w+ff ^ 7+7 



AT 

v 1] rf- {1 - 1/n)H * dM 1] r?. (3.56) 



x ^[i-Cl-i/^-H/o^Mo-^oTl x TT r n[l-(l-l/n)H e d/p e ] 



In particular, if n = • • ■ = rjy = r < 1, we combine (3.43) and (3.56) to obtain 
E[(L(x,T)-L(y,T)) n ] 

<c n _ \x ~ y\ ,n (nl)^ T e=i Hed/pe+He ° 1+1 , r "( ,v -( 1 - 1 / ,l )ELi ff " i / K "- H 'i) 7 ' 

^ < 24 («!) A '~ /3T+(1+ffT)7 N " y\ n ~<r n{ ^- H ^\ (3.57) 

The last inequality follows from the fact that Hi < H T and Lemma 3.4. This finishes the proof of Lemma 
3.10. □ 



Now we are ready to prove Theorem 3.1. It is similar to the proof of Theorem 4.1 in [33] and we include 
it for the sake of completeness. 
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Proof of Theorem 3.1. Let I E Abe fixed. For simplicity, we will assume / = [77, 1]^ for some 77 > 0, say, 
Tj = 2e (cf. (2.6)). It follows from Lemma 3.10 and the multiparameter version of Kolmogorov's continuity 
theorem (cf. [17]) that, for every T G A such that T C I, B H has almost surely a local time L(x,T) that is 
continuous for all x G M. d . 

To prove the joint continuity, observe that for all x, y G W 1 and s,t G /, we have 

E[(L(x, [r,, s}) - L(y, [n, t])) n ] < 2»- 1 {E[(L(x, [ V> s}) - L(x, [77, t])) n ] 

+ E[(L(x,[ V ,t])-L(y,[r ] ,t})) n }}- (3-58) 

Since the difference L(x, [rj, s]) — L(x, [17, t]) can be written as a sum of a finite number (only depends on N) 
of terms of the form L(x,Tj), where each Tj G A is a closed subinterval of / with at least one edge length 
< \s — t\, we can use Lemma 3.7 and Remark 3.9, to bound the first term in (3.58). On the other hand, 
the second term in (3.58) can be dealt with using Lemma 3.10 as above. Consequently, for some 7 £ (0, 1) 
small, the right-hand side of (3.58) is bounded by (|x — y\ + \s — t|)™ 7 , where n > 2 is an arbitrary even 
integer. Therefore the joint continuity of the local times follows again from the multiparameter version of 
Kolmogorov's continuity theorem. This finishes the proof of Theorem 3.1. □ 

We end this section with the following two technical lemmas, which will be useful in the next section. 

Lemma 3.11. Under the conditions of Lemma 3.7, there exist positive and finite constants c 3 26 and c 3 27 , 
depending on N, d, H and I only, such that the following hold: 

(i) For all a G I and hypercubes T = [a, a + (r)} C / with edge length r€ (0,1), 16 R d and all integers n>l, 

E[L(x + B H (a),T) n ]<cl 2e (nVj N -^r n ^, (3.59) 

where /3 t = N-t- H T d + £J =1 H T /H e . 

(ii) For all a G I and hypercubes T = [a, a + (r)] C I, x, y G M. d with \x — y\ < 1, all even integers n>l and 
all 'ye (0,lAf£), 

E[{L(x + B H {a) 1 T) - L(y + B H \a),T)) n ] < c n i 27 {n\) N -^ +{1+H ^\x - ( 3 , 60 ) 

Proof. For each fixed a € I, we define the Gaussian random field Y = {Y(t), t G M^} with values in M. d by 
Y(t) = B H (t) - B H (a). It follows from (3.1) that if B H has a local time L(x, S) on any Borcl set S\ then Y 
also has a local time L(x, S) on S and, moreover, L(x + B H (a),S) = L(x, S). With little modification, the 
proofs of Lemmas 3.7 and 3.10 go through for the Gaussian field Y. Hence we derive that both (3.59) and 
(3.60) hold. □ 

The following lemma is a consequence of Lemma 3.11 and Chebyshev's inequality. The proof is standard, 
hence omitted. 



Lemma 3.12. Under the conditions of Lemma 3.7, there exist positive constants c 3 28 ,c 3 29 ; b\ and bi > 
( depending on N , d, I and H only), such that for all a E I, T = [a, a + (r)] with r G (0, 1), x G M. d and u > 1 
large enough, we have 

P{L(x + B H (a),T) >c 3 28 r^u N ~^} < exp(-M) (3.61) 
and for x, y G M. d with \x — y\ < 1 and 7 G (0, 1 A 

P{|L(x + B H (a), T) - L(y + B H (a),T)\ 

> c a ,2o\x - y y r ^- H ^u N -^ +(1+H ^} < cxp(-M). (3.62) 
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4. Holder conditions for the local times 

In this section wc investigate the local and uniform asymptotic behavior of the local time L(x,T) at x and 
the maximum local time L*(T) = max x€R d L(x,T) as diam(T) — > 0. The results are then applied to study 
the sample path properties of B H . 

4-.1. Holder conditions for L(x,») 

By applying Lemma 3.12 (more precisely, (3.61) with a = 0) and the Borel Cantelli lemma, one can easily 
derive the following law of the iterated logarithm for the local time L(x,-): If (3.26) holds for some t£ 
{1, . . . , N}, then there exists a positive constant c 4 l such that for every x £ M. d and t £ (0, oo)^, 

L(x,U(t,r)) 

hmsup V / ' " <c 41 , 4.1 

where U(t,r) is the open ball centered at t with radius r and fi(r) = r@ T (\og\og(l/r)) N ~P T . It would be 
interesting to prove the lower bound in (4.1). For such a result for the local times of a one-parameter 
fractional Brownian motion, see [23] . 

It follows from Fubini's theorem that, with probability one, (4.1) holds for Ajv-almost all t£ (0,00)^. 
Now we prove a stronger version of this result, which is useful in determining the exact Hausdorff measure 
of the level set. 

Theorem 4.1. Assume that d < ^2t=i ~W~ e ■ Let t £ {1, . . . ,N} be the integer such that (3.26) holds and 
let I £ A be a fixed interval. For any fixed x £ M. d , let L(x,-) be the local time of B (t) at x which is a 
random measure supported on the level set (B H )^ 1 (x). Then there exists a positive and finite constant c i2 
independent of x such that with probability 1, 

L(x,U(t,r)) , in . 
limsup V V' " <c 4 2 4.2 

holds for L(x, ■) -almost all t£ I, where ip\ (r) = r^ T (loglog(l/r)) Ar_,3T . 

Proof. The method of our proof is similar to that of Proposition 4.1 in [30]. For every integer k > 0, we 
consider the random measure Lk(x,») on the Borel subsets C of I defined by 



exp(-|S-+i<e,£ H (i)-*>)d£dt. (4.3) 



2k 

Then, by the occupation density formula (3.1) and the continuity of the function y 1— ► L(y, C), one can verify 
that almost surely Lj~(x, C) — > L(x, C) as k — > 00 for every Borel set Cel. 

For every integer m > 1, denote f m (t) = L(x, U(t, 2~ m )). From the proof of Theorem 3.1 we can see that 
almost surely the functions f m (t) are continuous and bounded. Hence we have almost surely, for all integers 
m, n > 1, 

[f m (t)) n L(x,dt)= lim f[f m (t)] n L k (x,dt). (4.4) 
fc— *°° J 1 

It follows from (4.3), (4.4) and the proof of Proposition 3.1 of [25] that for every positive integer n > 1, 
E f[f m (t)] n L(x,dt) 



2n 



If I expf-ifj^,^)] 
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x Eexp^i^(u J ',.B H (a J ')}j duds, (4.5) 

where u= (it 1 , . . . , u n+1 ) € R( n+1 ) d and s = (t, a 1 , . . . , a"). Similar to the proof of (3.27) we have that the 
right-hand side of Eq. (4.5) is at most 

<c" 4 (n!) Af -' 3 -2-" m ^, (4.6) 



Uu(t,2-™ r ldetCov(B»(t),B*(si),...,B»(s-))y/2 



where c 4 4 is a positive finite constant depending on TV, d, H, and / only. 

Let 7 > be a constant whose value will be determined later. We consider the random set 

I m (w) = {tel: f m {t)>l^i{2- m )}- 

Denote by /i w the restriction of the random measure L(x, ■) on /, that is, ^(E) = L(x,E(lI) for every Borel 
set £cRf Now we take n — [logmj , where \_y\ denotes the integer part of y. Then by applying (4.6) and 
Stirling's formula, wc have 

[m(2 -m)]n 
c?>! )^ 2 -^ 

- 7 n 2 -mn/3 T (l g m )n(JV-/3 T ) - ' » ' 

provided 7 > is chosen large enough, say, 7 > c 4 4 e 2 := c 4 2 . This implies that 
E[ f^f, u (I m ) j <co. 

\m=l / 

Therefore, with probability 1 for /i w almost all t £ I, we have 

L(a:,E/(i,2- m )) „, 
Innsup ' i <c ii2 . (4.8) 

Finally, for any r > small enough, there exists an integer m such that 2~ m <r< 2~ m+1 and (4.8) is 
applicable. Since <^i(r) is increasing near r = 0, (4.2) follows from (4.8) and a monotonicity argument. □ 

Theorem 4.2. Assume that Yl^—i 7)7, > d and I £ A. Then there exists a positive constant c 4 5 such that 
for every x € M. d , 

ip 1 -m{{B H y\x)r\I)>c /L>5 L(x,I) a.s., (4.9) 
where tp± -m denotes the tpi -Hausdorff measure. 

Proof. As in the proof of Theorem 4.1 in [30], (4.9) follows from Theorem 4.1 and the upper density theorem 
of [26]. We omit the details. □ 



4-2. Holder conditions for £*(•) 



The following theorem establishes sharp Holder conditions for the maximum local times L*{T) = 
sup^ggd L(x, T) of fractional Brownian sheets as diam(T) — > 0. Similar results for Brownian motion and 
some other random fields have been obtained by several authors. See, for example, [12, 16, 20, 30]. 
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Theorem 4.3. Let B H = {B H (t),t £ R^} be a fractional Brownian sheet in M. d with index H = 
(Hi, . . ., Hn). We assume that there exists t £ {1, . . . , N} such that H\ = • • • = H T and Hid < r. Then, 
there exist positive constants c 4 6 and c 4 7 such that for every s £ /, 

lim sup ^* ( J S ~ (r) ' S + ^ < c 4 „ a.s. (4.10) 



and 



L*([s-{r),s+(r}]) 

limsupsup „ ,„ = . „ , < c , 7 o.s. (4.11) 

r ^o a6 j ^"^(logr-- 1 )^ ~ 4,7 



For proving Theorem 4.3, we will make use of the following lemma, which is a consequence of Lemma 2.1 
in [28] and Lemma 8 in [3]. 

Lemma 4.4. Let B H = {B (t), t £ R^} be a fractional Brownian sheet in R d with index H = (Hi, . . ., H^) 
and let I £ A be fixed. Then there exist positive constants c 4S and c 4 9 such that for all s £ I , T = [s, s + (h)] 
with h £ (0, 1) and all u> c 4 8 h Hl , we have 

p{sup|£ ff (i) - B H (s)\ > < cxp (- - U K2H } j ■ (4.12) 

Proof of Theorem 4.3. As in [12, 20, 30], the proof of Theorem 4.3 is based on Lemma 3.12 and a chaining 
argument. Hence we will only sketch a proof of (4.10), indicating the necessary modifications. 

Let g(r) = r N ~ Hld (\og\ogr~ 1 ) Hld for r > small enough. In order to prove (4.10) it is sufficient to show 
that for every s £ /, 

limsup ^r-%- <c 410 a.s., (4.13) 

where C n = [s, s + (2~ n )] for n > 1. 

We divide the proof of (4.13) into four steps. 

(a) Pick u = 2~ nHl ^/2c 4 9 logn in Lemma 4.4, we have 



»{ sup \B H (t) - B H (s)\ > 2-" Hl ^/2c 4 9 logn} < exp(-21ogn) = n~ 2 . (4.14) 



l teC„ 

Hence the Borel-Cantelli lemma implies that a.s. 3n% = ni(uj) such that 



tec„ 



sup\B H {t)-B H {s)\<2~ nHl J2c i Aogn, foralln>ni. (4.15) 



(b) Let 6 n = 2-" ffl (loglog2")-( 1+Hl ) for all n> 1, and define 



G n = {x £ R d : |x| < 2~ Tlffl v/2c 4jB logn with x = 9 n p for some p £ Z d }. 
Then, at least when n is large enough, the cardinality of G n satisfies 

%G n <c 4A1 (\ogn)( 2+H ^ d . (4.16) 
It follows from (3.61) that for any constant c > and integer n large enough, 

p\maxL{x + B H (s),C n )>c Hld g{2~ n )} <c 41 .,(\ognY 2+Hl)d n- cbl . (4.17) 
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(Note that (3 T = N — Hid under the assumptions of Theorem 4.3.) By choosing c= 2&7 1 in (4.17) we see 
that the right-hand side of (4.17) is summable. Hence, the Borel-Cantelli lemma implies that almost surely 
3n2 = n2(u) such that 

maji L(x + B H (s),C„) < (2b7 1 ) Hld g(2- n ), foralln>n 2 . (4.18) 
(c) Given integers n, k > 1 and x £ G n , we define 

F(n, k, x) = \ y £ R d : y = x + 



k \ 

^ £j -2-^ £j £{0,l} d for l<j<k\. 

j=l J 



A pair of points yi,?/2 £ F(n, k,x) is said to be linked if yi — yi = n e2 k for some e £ {0, l} d . We choose 
7 > small such that (3.62) in Lemma 3.12 holds, and then choose S > such that S(Hid+ (1 + #1)7) < 7. 
Consider the event B n defined by 

00 

B «= U U \J{\L(yi+B H (s),C n )-L(y 2 + B H (s),C n )\ 

i6G„ k=l J/l,I/2 

>2-^ N - H ^- H ^\ yi -y 2 p(c2 5k \ogn) Hld+{1+Hlh }, (4.19) 

where (J y9 signifies the union over all the linked pairs 2/1,2/2, and where c > is a constant whose value 
will be chosen later. 

Note that H T = Hi, by (3.62) we derive that for n large enough, 

00 

P{B n } < c 4 1 3(logn) (2+Hl)d ^2( d+1 ) fc exp(-cfe 2 2 ifc logn) 

fe=i 

<c il4 (logn)^ +H ^ d n~ cb \ (4.20) 

In the above the last inequality follows from the fact 

00 

2 (d+1)fc exp(-x2 5k ) <e- x , Vx > large enough. 

fe=i 

Hence, by choosing c= 2&2 -1 in (4.19), the Borel-Cantelli lemma implies that almost surely, B n occurs only 
finitely many times. 

(d) Fix an integer n together with some y € M. d that satisfies \y\ < 2~ nHl ^/2c 4 9 logn, we can represent y 
in the form y = lim^oo y^ with 

k 

y k = x + 6 n J2zj2~ j > ( 4 - 21 ) 
3=1 

where yo = x £ G n and £j £ {0, l} d for j = 1, . . . , k. 

Since the local time L is jointly continuous, by expansion (4.21) and the triangular inequality, we see that 
on the event , 

\L(y + B H {s),C n ) - L(x + B H (s), C n )\ 

OO 

< J2\L(Vk + B H (s),C n ) - L(y k .i + B H (s),C n )\ 

k=l 
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oo 

<J22'^ N - H ^- H ^\y k -y k ^(2b^ 1 2 Sk logn) Hld+{1+Hl} '' 
fc=l 

<c 4 , 15 p(2-"). (4.22) 
We combine (4.18) and (4.22) to get that for n large enough, 

sup L(x + B H (s), C n ) < c 4 , le5 (2-"). (4.23) 

\x\<2-" H i v /2c ig logn 

That is 

sup L(x,C n )<c il6 g(2- n ). (4.24) 

\x-BH( S )\<2-" H i^/2c i glogn 



Since L*(C n ) = sup{L(;r, C„): a; G 5 ff (C n )}, (4.13) follows from (4.24). This proves Theorem 4.3. □ 

The Holder conditions for the local times of fractional Brownian sheets are closely related to the irreg- 
ularity of the sample path of B H \t). To end this paper, we apply Theorem 4.3 to derive results about the 
degree of oscillation of the sample paths of B H (t) 1 which greatly improves Theorem 3 of [3]. 

Theorem 4.5. Let B H = {B H (t),t G R+} be an (N,d) -fractional Brownian sheet and let I E A be a fixed 
interval. Then there exists a constant c 4 17 > such that for every s G /, 

l™inf sup _L B S ~ f-nlg, ^ <Vr a*. (4.25) 



t£U(s,r) 



r) ^(loglogr" 1 )- 



and 



\B H (t) — B H (s)\ 

liminfinf sup ' }' -, N „ > c A a.s. (4.26) 

tGi tec/(s>r) ^(logr- 1 )-^ " 4 ' 17 



/v 



In particular, the sample function B (t) is almost surely nowhere differ entiable in (0,oo) 

Proof. It is sufficient to prove the results for d = 1. Note that Hi < 1, Theorem 4.3 is always applicable for 
d = 1 with r = 1 . For any interval Q G A, we have 



A 



(Q)= I L{x,Q)dx<L*(Q) x sup |flf(u)-B^(«)|. (4.27) 

JB»{Q) u.veQ 



By taking Q = ^(s^) we see that (4.25) follows immediately from (4.27) and (4.10). Similarly (4.26) follows 
from (4.27) and (4.11). □ 
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